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^ ■ Abstract 
0\ 

A microscopic treatment of cumulative phenomena based on perturbative QCD 
^ ' calculations of the corresponding quark diagrams near the threshold is presented. 

To sum all diagrams like these the special technique based on the recurrence rela- 



tions was developed. The x-behaviour of the nuclear structure function F2{x) in the 
cumulative region x > 1 was found to be roughly exponential, governed by an ef- 
fective coupling constant, which depends on the QCD coupling constant and quark 
, mass. Two mechanisms for cumulative particle production, direct and spectator 

ones, were analysed. It was shown that due to final state interactions the leading 
Ps| ' terms of the direct mechanism contribution are cancelled and the spectator mech- 

. anism is the dominant one. It leads to a smaller slope of the particle production 

\^ I rates compared to the slope of the nuclear structure function in the cumulative re- 

0^ ' gion, in agreement with the recent experimental data. The slope difference is due to 

additional multiple interactions between nuclear and projectile partons which enter 
Qh! the spectator mechanism for the cumulative production. The different versions of 

hadronization mechanisms of the produced cumulative quarks into hadrons are also 
5^ ■ discussed. 

> 

X. 

S ; 1 Introduction 

One of the most interesting fragmentation processes is the production of fast hadrons 
in hA interactions in the nuclear fragmentation region with the value of the Feynman 
scaling variable: x > 1, i.e. in the region in which the production of this fast hadron is 
kinematically prohibited for the hN interaction. Here x is the longitudinal momentum 
per nucleon of the fast moving nucleus carried away by the produced particle. It means 
that a produced particle carries a momentum much greater than the average one of a 
constituent nucleon. The nuclear fragmentation processes like these usually are referred 
as the cumulative ones. For the deep inelastic I A scattering the cumulative region is x > 1 
in terms of Bjorken scaling variable defined for individual IN interactions. 

The theoretical investigations of the nuclear fragmentation processes in this cumulative 
region are of great physical importance as they enable to study multi-nucleon short-range 
correlations in nuclei (so called fluctons) and to get the information on the high-density 
nuclear matter clusters which always are being in nucleus. 
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From the modem point of view these compact multi-quark dense hadronic matter 
clusters in nuclei can be also considered as the quark-gluon plasma clusters at zero tem- 
perature which investigation is of interest also in the light of recent efforts to detect the 
indications of high temperature quark-gluon plasma formation in heavy-ion collisions on 
colliders at super high energies. 

An adequate treatment of the cumulative effect requires the quark language |^ . A 
phenomenological treatment of cumulative phenomena based on nuclear quark distribu- 
tions and QCD evolution equations was proposed in It uses an ad hoc relation between 
the structure functions and particle production rate, based on the nucleonic (or fluctonic) 
picture of the nucleus rather than the quark one. Starting from the quark picture of 
hadrons the famous quark counting rules have been proposed which, in particular, 

describe the threshold behaviour of the structure functions. In it was noted that these 
rules also apply to nuclei at the deep cumulative threshold x ^ A. In [|, Q it was stressed 
that the cumulative phenomena involve a particular type of contributions, the so-called 
intrinsic diagrams, in which several partons of the nucleus cooperate to determine the 
intrinsic hardness of the nuclear wave function. 

However the quark counting rules alone do not allow to find the cumulative probabil- 
ities. To calculate these, one needs besides to know the coefficients of the quark counting 
rules behaviour, which determine the relative weights with which different contributions 
combine to give the total probability, as well as the overall normalization. Recently we 
have proposed 0, H a microscopic model for the theoretical description of the cumulative 
phenomena based on perturbative QCD calculations of the corresponding quark-parton 
Feynman diagrams near the thresholds and calculated these coefficients in the framework 
of the proposed model. By means of the specially developed technique based on the recur- 
rence relations all diagrams of that kind were summed. As a result, we obtain formulas 
for the nuclear structure function and particle production rate at x > 1, which allow 
direct comparison with experiment. The proposed method of the summation of intrinsic 
diagrams at cumulative thresholds based on the recurrence relation technique are also 
discussed in the present paper. 

2 General formulae 

2.1 The nuclear structure function at x > 1 

Consider the virtual 7-nucleus scattering cross-section in the system, where the nucleus 
is moving fast along the z-axis with the longitudinal momentum per nucleon pz- As a 
reasonable first approximation, we may treat the nucleus as a colllection of = 3A 
valence quarks, which, on the average, carry each momentum xoPz with xq = 1/3. In 
the following we assume that the quark longitudinal momentum distribution is sharp 
enough, so that in all places, except in the quark wave function itself, we can safely put 
the initial quark longitudinal momentum equal to xoPz- The probing 7 quantum will 
see a "active" quark with its momentum much greater than xoPz only if this quark has 
interacted with several other quarks ("donors") and has taken some of their longitudinal 
momenta. The maximum possible value of the longitudinal momentum thus accumulated 
is NxqPz = Apz- It corresponds to A^ — 1 interactions with all other quarks, whose 
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Figure 1: The deep inelastic 'jA scattering 

longitudinal momentum has become equal to zero each. It is well-known that interactions 
which make the longitudinal momentum of one of the quark equal to zero may be treated 
by perturbation theory [Q. These considerations were used in numerous calculations of 
hadronic structure functions near a; = 1, which revealed that the results depend heavily 
on the quark distribution in the initial hadron (see e.g [Q). In the case of a nucleus the 
calculations in principle remain the same. With a number of quarks considerably larger, 
their distribution in x evidently becomes much narrower, so that we may hope that the 
results would be less sensitive to its form nor to its spin-colour structure. On the other 
hand the number of relevant diagrams rises tremendously. Their summation presents the 
main techical difficulty in the nuclear case. 

In general form the diagram for the deep inelastic scattering is shown in Fig. |l[ 
In the following we use the light-cone variables p± = {po ± pz) / V^- The momentum of 
the active quark will be denoted as k and its scaling variable is ,^ = k^/p+ ~ kz/pz- We 
shall study a general case when n — 1 quarks aquire their momenta equal to zero as a 
result of interquark interactions, giving their momentum to the active quark. Evidently 
with n — 1 momentum exchanges the ^ maximal value is uxq. The upper blob in Fig. |I] 
represents the active quark's structure function. For simplicity we shall treat quarks as 
scalar particles. Then the contribution of the diagram of Fig. || to the nuclear structure 
function ^2"^'' is given by 



[x,q 



{di/0m)f2{x/i,q' 



where D{^)/^ has the meaning of the probability to find in the nucleus a quark with the 
longitudinal momentum kz = ^Pz, /2 is the quark structure function and x = q^/2qp. In 
the cumulative region 1 < x < ^ < A. 



2.2 Cumulative particle production 

We now turn to the production of cumulative particles which have x > 1 in the system 
where the nucleus is moving fast along the z- axis. All the contributions can be divided 
into a direct part, in which the projectile interacts with the created cumulative quark 
(Fig. ^) and which is a straightforward generalization of Fig. 0, and a spectator part in 
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Figure 2: Cumulative particle production, the direct mechanism. 



which the projectile interacts with other quarks. The direct contribution is evidently given 
by the formula analogous to (|ID. We shall limit ourselves with the inclusive cross-section 
integrated over the transverse momentum Ia{x) = xda/dx. The direct part is then 

iTix) = I {di/ODim^/O (2) 

J X 

where again D{^)/^ is the probability defined in the preceding subsection to find in the 
nucleus a cumulative quark with scaling variable ^, h{x/^) is the inclusive cross-section 
with this quark as a target, x = k+/p+ ^ kz/pz and ^ = k'_^_/p+ ~ K/Pz- Again in the 
cumulative region 1 < x < ^ < A. 

The spectator part has a much more complicated structure. We have analysed its 
contribution in details in 0, We have found that near the thresholds significant 
cancellations occur in this part, which lead to the domination of the diagrams, with each 
donor quark interacting with the projectile . As a result the contribution of the spectator 
mechanism has not such a simple convolutional form (0) as in the direct case. With the 
quasi-eikonal parametrization of the partonic amplitude chosen to account for diffraction, 
we found ||^, ^ for the contribution of the spectator mechanism to the cumulative particle 
production in the case when n — 1 donors transfer their momenta to the active quark: 

ly^ix) = Dn{x)Cp' J d'B[47rm'j{B)r-' (3) 

with 

JiB) = CpUpJ dVb'X\b)pp{b')\a{B + b'~b)- a{B + 6') I' (4) 

where 

X{b) = Ko{mb)/2Tr (5) 

Here up is the mean number of partons in the projectile, pp is the one parton distribution 
normalized to unity and Cp is the quasi-eikonal factor for the projectile. For example, for 
NN-interaction Cfj- = 1 + ct'^nn/^nn- The a is the partonic amplitude of the interaction 
between a parton of the projectile and a donor quark from the nucleus. 
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Figure 3: Formation of the hard parton component of the flucton in the nucleus. 
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Figure 4: An example of the diagram which contribute to Bn- 

3 The cumulative quark production probability 

We see that both the nuclear structure function (P and the particle production rate (||J^) 
in cumulative region are determined by -D(^)/^ - function which has the meaning of the 
probability to find in the nucleus a quark with the scaling variable ^ > 1. 

Study the general case when n — 1 quarks aquire their momenta equal to zero as a 
result of interquark interactions, giving their momentum to the active quark (Fig. 
The scaling variables for quarks are defined as = According to our assumption 

initially in the nucleus all Xi ~ xq and the distribution in Xi is sharp. The momentum of 
the active quark is denoted as = fci and its scaling variable is = xi. Evidently with 
n — 1 momentum exchanges its maximal value is nxo and when xi = C, nxo then all 
Xi ^ for i = 2,3, ...,n. Particulaly on this fact the use of the perturbation theory is 
based §. 

The block Bn describes the formation of the hard parton component of the flucton as 
a result of n — 1 interquark interactions. Due to the quark counting rules Q one has 
DniO ~ {fT'^o — ^)^"~^ when ^ — ^ hxq. But as we have mentioned in the introduction for 
to calculate the cumulative quark production probability one needs besides to know the 
coefficients of the quark counting rules behaviour, which determine the relative weights 
with which contributions with different n combine to give the total probability, as well 
as the overall normalization. So unlike we are trying to calculate theoretically the 
structure function of the flucton at least its hard component. 

In the perturbative approach we find Feynman diagrams for the i?„ of the type shown 
in Fig. H. The shown diagram is only an example. It contains a particular set of interquark 
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Figure 5: Decomposition of the Feynman diagram on the time ordered diagrams. 

interactions, in which donors successively give their momentum to the active quark and 
become soft. However other type of interquark interactions are also possible: say, the 
n — 1 donors may give their momenta successively to each other, the last giving all the 
accumulated momentum to the active quark (see second diagram in Fig. ^). On these 
diagrams dash vertical lines refer to one gluon exchange. One can choose the Coulomb 
gauge in which the transversal part is damped at low Xi and the dominating Coulomb 
part depends only on the scaling variables of the participant quarks before (xi,X2) and 
after {x^^x^) the interaction (see Fig. |^). If we neglect the colour then 

V = Ana^xi + x[){x2 + X2)/{xi — x[Y = A-koV (6) 

with a the interaction constant. To approximately take the colour into account we average 
over the quark colour variables, which introduces a factor -\/2/3 into (||). 

We are using the time ordered pertubation theory in the light cone variables. In the 
energy denominators of the intermediate states the contributions from soft final quarks 
dominate as their "energeis" Ei = {m? + k^_|_)/xi rise with Xj — for i = 2, 3, n. Thus 
the contribution of the diagram in Fig. ^ will be proportional to [E2{E2 + E^)...{E2 + ... + 
En)Y^ or after symmetrization on X2, to [E2E^...Ej^~'^ . 

The advantage of the time ordered pertubation theory is that the energy denomi- 
nators have the same form for any diagram in the threshold limit and coincied after 
symmetrization on a;2,...,a;„. To illustrate, we have presented in Fig. ^ one Feynman 
diagram as the sum of two time ordered diagrams. The Feynman diagram is propor- 
tional to [£'2-^4 (-^2 + -E3 + -^4)]"^ and the time ordered diagrams are proportional to 
[^2(^2 + Ei){E2 + Ei + E^)]-^ and to [^4(^4 + ^2)(^4 + ^2 + ^3)]"^ respectively in 
the ^ — >■ 4xo limit. After symmetrization on X2,...,x„ the contributions of the last two 
diagrams both become proportional to [E2E^Ei\~'^ . Moreover in this case they are equal 
which leads to the possibility to take into account only one of them multiplied by the 
"time ordering" factor 2. 

As result the -B„ can be presented in the form 

" " {xoY-^ E2E,...E^ 

The denominator leads after integration over X2, x„ to the quark counting rule behavior 
for the cumulative quark production rate from the n-quark flucton. But our aim is to 
calculate the coefficient in this behavior, which determines the relative weights with which 
contributions on different thresholds combine to give the total probability. For to do 
this we have to calculate the sum of all possible Feynman diagrams of the type shown 
in Fig. ^. In (^ the X„ denotes this sum after the extraction of the common factor 
(47ra)"-i(a;o)2~"[S2^3...^n] 
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Figure 6: The diagrams which contribute to X„ for n = 6 as an example. All are in 
units Xq. 
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After this extraction one can go to the limit Xi = C, = uxq and Xj = for i = 2, n 
in the rest of the diagram, which leads to the diagram rules in Fig. |^ for Xn- All Xi in 
Fig. are in units xq. We have also to multiply each diagram on the proper "time 
ordering" factor for to take into account the different time ordering possibilities for gluon 
exchanges. 

The main problem here is that the number of these diagrams grows rapidly with n, 
as n\. For to sum all these diagrams we derive a recurrency relation. The derivation is 
based on the important observation that only diagrams with all intermediate x greater 
than zero should be retained. Separating the last exchange, which raises the momentum 
of the active quark up to its final value we present X„ as a sum over k = l,2,...,n — 1 
of the diagrams shown in Fig. ||, where it has been assumed that Xi = 1 by definition. 
Times of /c — 1 gluon exchanges in the upper blob X^. relative to those in the lower blob 
Xn-k may be arbitrary, which leads to a time ordering factor C^Zl- Taking into account 
the diagram rules in Fig. |^ and denoting /„ = X„/n! we find that Fig. ^ translates into 
a recurrency relation for /„: 
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with the initial condition /i = 1. The recurrency relation enables easy calculate /„ 
for an arbitrary n starting from /i = 1. For large n (El) evidently admits asymptotical 
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Figure 8: The diagrammatic recurrency relation for X„. 
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solutions of the form 

/„ ~ [(6/5)n + o{n)] exp(— an) (9) 
where a is arbitrary. Numerical studies reveal that with /i = 1 

a = 0.24421... 

and also show that the asymptotical expression appproximates the true solution quite 
well starting from n = 3, i.e. for all physically interesting values. 

So we have calculated the X„ = /„«! and consequently the Bn which describes the 
formation of the hard parton component of the flucton. Now from diagram in Fig. |^ we 
can calculate the D{C,)/C, - the probability to find in the nucleus a hard cumulative quark 
with scaling variable C,- 

The block i?„ is not dependent on the momentum variables of the initial quarks. As 
the consequence we find that all donor quarks have their coordinates relative to the active 
one equal to zero. 

Now assume that on the average quarks are homogeneously distributed inside nucleons 
(three in each of them) . In this picture the probability w^^^ to find n — 1 quarks within the 
distance Ri from the active one, where Ri is the nucleon radius, is equal to the probability 
to find within the same distance n nucleons, where n = [(n — l)/3] + 1 (i.e the entire part 
of {n — l)/3 plus one nucleons): 

wi^^ = 4""^ (10) 

If the nucleons, in their turn, are homogeneously distributed inside the nucleus of the 
finite volume Va = AVq then the latter probability will be given by 

4""^ = crMv,/VAr-' (11) 

where Vi is the nucleon's volume. Actually, as well-known, this approximation is rather 
poor due to the rise of the nuclear density towerds its center. To account for this fact we 
introduce a factor g„ to be calculated from a realistic nuclear distribution. In particular 
for the Saxon- Woods distribution 

p(r) = po/(l + exp(r - i?^)/ri)) (12) 

with Ra = A^^^ 1.17 fm the nuclear radius and ri = 0.51 fm, for A >> 1 

1 n— 1 

qn = l- {3n/RA) = 1 - 1.31 ^(1/A;) (13) 

k=l k=l 

An explicit calculation gives for the cumulative quark distribution 



A.(o = 3^c.r_;A-^5^j ^[^j (14) 

where n — 1 = [{n — l)/3] is the "cumulative number" (i.e. the number of nucleons 
in the flucton n minus one) and = {nxo — O/^o- AH nontrivial A dependence of 
the distribution D is contained in the nuclear factor qn- It rises with A approaching its 



8 



asymptotic value 1. The rise is quicker for higher n. Studies of the cumulative particle 
production rate in the nucleon picture, where the same A- dependence emerges [ p!0[| , 
revealed that it agrees rather well with the experimental observations. 
With d) we finally find 

DniO = ^ (^4l)y^- (^)"" - l)-^/^"-^^(A„)A--^ (15) 

where, 

7 = aV36mVi exp 2a (16) 

and we have assumed A » 1. The behavior of Dn{C,) near the threshold A„ = is 
determined by the factor A^"~'^, which corresponds to the well-known quark counting 
rules 0. The magnitude of the coefficient is determined mainly by the factor 'y^~^. As 
we shall see later, the values of 7 extracted from the known experimental data on nuclear 
structure functions in the cumulative region result quite small. Because of this the power 
factor 7""^ ensures that the coefficient in the ^ dependence of ([T5|) drops very fast with the 
growth of n. Down from the treshold the contribution Dn{C) rises as a power. This rise has 
evidently to stop at some distance from the threshold where our approximations become 
invalid and terms with higher powers of A„ become important. We expect, however, that 
the overall order of the contribution remains to be governed by the small factor 7""^ 
coming from n — 1 gluon exchanges. 

As a result, due to the smallness of 7, the contribution from diagrams with a given n 
will always dominate over those with greater n's, as soon as we slightly move down from 
the threshold point ^ = uxq. Therfore we may safely assume that Dn{C) given by ([T^) or 
(|l5| ) represents the true quark distribution for ^ starting from the threshold ^ = nxo down 
to the point ^ = {n — 1)xq, which is the threshold for the diagrams with n — 2 donors, i. 
e. in the region 

{n — l)xo < ^ < tixq, Xq = 1/3 (17) 

Alternativelely, to avoid discontinuities, we may take the total distribution D{C,) as a sum 
of Dn{C) over all possible n's with an additional assumption that A„ = 1 for ^ smaller than 
(n — 1)2:0, which means freezing the contribution D„ at the threshold for the contribution 
Dn-i- As a result we have the roughly exponential behaviour of D{^) in the cumulative 
region, governed by an effective coupling constant, which depends on the QCD coupling 
constant and quark mass. 

4 Numerical calculations, discussion and the com- 
parison to experimental data 

Both convolutional formulas (|l],0) lead to the same x dependence for the structure function 
F2{x) and for the direct mechanism production rate J^*''(x) in the region x > 1. It is 
roughly an exponential in x 

If'{x) ~ F^{x) ~ exp(-6oa;) (18) 
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where the slope bo is determined by the QCD couphng constant and quark mass. 

The spectator mechanism besides involves interactions between partons of the 
projectile and target nucleus. It was shown in that each donor quark has to interact 
with the projectile. As a result the spectator contribution I^i^) also behaves in x as an 
exponential but with a different slope 

r/{x) ~ exp(-6,x) (19) 

the slope bg depending not only on the QCD coupling constant and quark mass but also 
on the partonic amphtude. 

In IP we have found that final state interactions cancel the leading terms in the direct 
contribution, so that it becomes much less than the spectator one. Then particle produc- 
tion in the cumulative region indeed goes predominantly via the spectator mechanism. 
To be able to explain the experimental slope we have studied in P] the quasi-eikonal 
parametrization for the partonic amplitude. The spectator slope bg results very sensi- 
tive to the magnitude of the hadronic diffractive cross-sections. We have chosen the 
maximum possible value for the diffractive parameter for the case of the nucleon as a pro- 
jectile Cp = Cn = 1.5 (see (^)). Then the quasi-eikonal parametrization leads to values 
of the partonic amplitude a (^) considerably larger than without diffraction because of a 
stronger screening effect introduced by diffractive states. 

The results of the calculations are shown in Fig. 9 for the cumulative charged pion 
production on ^^^Ta, respectively. For comparison the corresponding nuclear structure 
function F.^{x) at a; > 1, is also shown, as well as the available experimental data from 

One clearly observes that the spectator mechanism, with a quasi-eikonal parametriza- 
tion of the partonic amplitude chosen to account for diffraction, leads to a considerably 
smaller slope of the production spectra (bg ~ 7^9) (|19D compared to the slope of the struc- 
ture function (|l^) in the region x > 1 (6o ~ 16), in a good agreement with experimental 
data. 

Thus our model correctly predicts two differents exponential in the cumulative particle 
production and in the nuclear structure function at x > 1. The difference is due to 
additional multiple interactions between projectile and target which enter the spectator 
mechanism for the cumulative production. 

A more phenomenological attempt to explain smaller slopes for the cumulative particle 
production is made in WM, where the existence of multiquark clusters in the nucleus and 



their properties are postulated and the quark-gluon string model formalism of ||T^ is used 
to calculate production rates. 

Note that at a nucleon level of analysis a similar approach to a description of cumula- 



tive phenomena was suggested in our paper ||16[ in which to find the cumulative particle 



production rates we used the fragmentation functions of multi-nucleon fluctons calculated 
in the framework of Regge approach. In a sense, the calculations of the probabilities "to 
slow" the quark system with some definite quantum numbers usually discussed in the 
quark-gluon string model approach [l^, |l8l is close to the idea of counting of the 



quantum numbers transferring from the flucton fragmentation region to the central one 
in triple reggeon like approach [ff6l. In the last case these quantum numbers transferring 
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in t-chanel determine the number and the type of Regge trajectories which have to be 
included in the "shoulders" of the triple like reggeon diagram. 

To calculate the hadron production rates in the cumulative region we need also to ana- 
lyze the process of fragmentation of one (or several) fast quark(s) into hadrons of different 
flavors. In principle the different versions of hadronization mechanisms of the produced 
fast cumulative quarks into colorless states (hadrons) can be suggested. Apparently that 
in the case of the production of cumulative protons the hadronization through the coales- 
cence of three cumulative quarks is favorable than the usual hadronization through one 
cumulative quark fragmentation into proton. This leads to the higher production rates 
for the protons compared with the pions but approximately with the same slopes of the 
protons and mesons spectra in agreement with experiment. 

The idea that in QCD a quark can hadronize by coalescing with a comoving spectator 
parton was suggested in the paper [|19|. It was used later for the description of the 
fragmentation of protons and pions into charm and beauty hadrons at large x pO| , It 
was shown that the coalescence or recombination of one or both intrinsic charm quarks 
with spectator valence quarks of the Fock state leads in a natural way to leading charm 
and beauty production. 

The conception of quark fusion for the description of hadronic fragmentation into vr^ 
and mesons was also used in the papers |23| . It was shown that the valence quarks 
of the projectile hadron play a dominating role in meson production in the projectile 
fragmentation region and while the x-distribution for vr"'", vr" and are the convolutions 
of the Uv, dv and Uy valence quarks with those of the corresponding antiseaquarks dg, Ug 
and Sg, the ^-distribution for K~ is that of the a sea quark Sg and an antiseaquark dg. The 
spin properties of these fragmentation processes have been also studied within framework 
of this approach p4[-|^. In particular it has been shown that the observed left-right 
asymmetry can be readily described in the framework of a relativistic quark model in 
which the observed vr"*" and 7i~ are respectively the fusion products of the valence quarks 
Uy and dy of p(t) and antiseaquarks dg and Ug. 

Note that both the intrinsic mechanism of the cumulative quark production when the 
quarks of several nucleons concentrated in one nuclear fiucton transfer their longitudinal 
momenta to the distinguished quark and the hadronization through the cumulative quarks 
coalescence break the QCD factorization theorem. 

This work is supported by the Russian Foundation for Fundamental Research, Grant 
No. 94-02-06024-a. 
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Figure 9: I a = is the calculated inclusive cross-section (per nucleon) for cumu- 

lative charged pion production on ^'^^Ta at y/s = 23.5 GeV (solid curve) and 1800 GeV 
(dashed one). A - the experimental data ||12| on the cumulative charged pion production 
on ^^^Ta by 400 GeV incident proton beam. F,^/A is the calculated nuclear structure 
function for the ^^^Ta at = 50 GeV^ (dashed curve) and at = 100 GeV^ 
(solid curve). □ and x - the experimental data |13[ on the ^^C structure function at 
g2 = 61 GeV^ and 150 GeV^, respectively. 
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